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Abstract 

We present a model independent procedure for extracting deep-inelastic structure function of 
"free" neutron from the electron - deuteron scattering with protons produced in the target frag- 
mentation region of the reaction. This procedure is based on the extrapolation of t, which 
describes the invariant momentum transfered to the proton, to the unphysical region corre- 
sponding to the mass of the struck neutron. We demonstrate that the impulse approximation 
diagram of the reaction has a pole at this limit with a residue being proportional to the "free" 
neutron structure function. The method is analogous to that of Chew and Low for extrac- 
tion of the "free" pion-pion and neutron-neutron cross sections from p(-ir,p)X and d(n, n)pn 
reactions respectively. We demonstrate that in the extrapolation the final state interaction am- 
plitudes are smooth functions of t and have negligible contribution in the extracted "free" 
nucleon structure function. We also estimate the range of the recoil nucleon momenta which 
could be used for successful extrapolation procedure. 



1 Introduction 

In spite of three decades of studies of the partonic structure of nucleons one still lacks a satisfactory 
knowledge of the relative d and u quark densities at large Bjorken x region. 

To extract these quantities two major approaches have been considered to date: one is the 
neutrino scattering off the proton at large x which allows one to probe separately the u and d distri- 
butions. The second is the extraction of parton distributions from both protons and neutrons using 
inclusive scattering from the hydrogen and deuteron targets. (In the future it would be possible 
also to use W ± production at LHC[ 1].) 
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While deep inelastic neutrino-proton scatterings lack adequate statistics, the inclusive electron- 
deuteron measurements suffer from significant nuclear effects an estimation of which involves the 
consideration of specific models for the Fermi motion and the EMC effect 121111 00. 

In our previous works O [SJ we outlined an alternative, model independent approach for 
extraction of "free" neutron deep inelastic structure function using semi-inclusive tagged neutron 
reactions, d(e, e', p)X, in which slow protons are detected in the target fragmentation region of 
the reaction. These considerations were incorporated in the experimental proposal which was 
approved recently at Jefferson Lab (JLab) @ l . Preparations for the experiment are currently 
under way and hence it is timely to elaborate our approach quantitatively. 

In our consideration of semi-inclusive deep inelastic scattering off the deuteron we focus in 
the region of relatively large x > 0.3 were coherence length is small and hence nuclear shadowing 
mechanism of the distortion of the nucleon spectrum discussed in [ 10 1 is negligible. 

The extraction procedure is based on the observation that due to a weak binding in the deuteron 
the singularity of the amplitude in the t = (pd—p P ) 2 -channel is much closer to the physical region 
of on-shell neutron than to all other singularities (the closest of which would be the pion production 
threshold). Hence one can in principle to continue analytically the scattering amplitude in t and 
find the residue at the pole of the struck neutron propagator at t — m 2 n . This is analogous to the 
Chew-Low procedure ifTTl for the extraction of the pion-pion scattering amplitude from the pion- 
nucleon data. It is worth emphasizing that in our case the analytic continuation is simpler since the 
elementary amplitude does not contain factors which go through zero at t close to the pole of the 
amplitude ( t = as compared to t — rnfy. Since procedures of extrapolation work well for the 
pion case, we expect them to be even more effective for the case of the tagged nucleon processes. 

The paper is organized as follows: In section II we outline the general framework of semi- 
inclusive deep inelastic scattering off nuclear targets, concentrating on the general properties of 
the scattering amplitude in the impulse approximation as well as beyond the impulse approxima- 
tion. In section III we study the analytic properties of the impulse approximation and the final 
state interaction amplitudes at kinematics of recoil proton extrapolated to the pole values of the 
struck neutron propagator. First we study the singular behavior of the impulse approximation am- 
plitude in the limit t — > m 2 n . Then we prove the loop theorem, which states that any additional 
loop in the scattering amplitude removes the singularity associated with the on-shellness of the 
struck neutron. In section IV we use a specific model to calculate the corrections to the impulse 
approximation due to the final state interactions. We observe that the deviations from the im- 
pulse approximation due to such reinteractions is mostly concentrated at recoil angles < 100-120°. 
Based on the considered final state interaction models we elaborate the extrapolation procedure 
and demonstrate that extracted "free" neutron DIS structure function is insensitive to the strength 
of the rescattering amplitude. This explicitly demonstrates the model independence of the extrap- 
olation procedure. We find that these procedure can be performed reliably if the tagged nucleons 
with momenta 50 4- 150 MeV/c are detected. 

'it is certainly of interest to study the reaction with tagged neutron as well which would allow a direct comparison 
of the scattering off a free and bound proton and hence identify the processes beyond the impulse approximation. 
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2 General Framework 



We consider semiinclusive deep inelastic scattering (DIS) off the deuteron: 

e + d-^e' + N + X 



(1) 



where nucleon N is detected in the target fragmentation region. We define E e as the initial energy 
of the electron and E' e , 6 e as the energy and scattered angle of the final electron, q = (z/, q) is the 
four momentum of the virtual photon, with v = E e — E' e and Q 2 = —q 2 . The recoil nucleon is 
described by four-momentum p s = (E s , p s ). We identify the masses of recoil and struck nucleons 
by m s and respectively, t = (pd — p s ) 2 determines the invariant momentum transferred to the 
recoil nucleon, where pd is the four momentum of the deuteron. 

The processes of Eq.([l]) with recoil proton can be used to extract the DIS structure function 
of the neutron. Since in this case the neutron is bound, it requires a careful treatment of off-shell 
effects in maximally model independent way. The approach we will discuss in this work is similar 
to one of Chew and Low ifTTIl who studied the issues of extracting n + + n° and n + n — > n + n 
cross sections from n + +p — > p + X and n + d — > p + n + n reactions respectively. In their analysis 
they observed that the analytical structure of the impulse approximation amplitude ( as in Eq.©) 
is such that it has a pole in nonphysical region of t corresponding to the one-mass-shell kinematics 
of the bound particles involved in the interaction. 

In the similar way we will analyze analytic properties of the scattering amplitude of the reac- 
tion (HJ) in the impulse approximation focusing on the issues related to the extraction of the on-shell 
DIS structure function of the neutron. 

First, we summarize the general formulae for the cross section of process ([T]): 



da 



Anat 



dxdQ 2 d 3 p s /E s xQ 4 



2 2 2 

x y m N . 



Q 2 



f t + + tan 2 (-))^cos((j))F^ L + cos (2(f)) F^ T 



Q 



.6 



D 



(2) 



where four independent nuclear structure functions, Ff T TL TT depend on Q 2 , x, a s , p st , with a s 



being light cone momentum fraction of the deuteron carried out by recoil nucleon. The 
latter is normalized in such way that the sum of a s + a = 2, where a is the similar quantity for the 
interacting nucleon. The Bjorken x = 2 ® and y = The z axis aligned in the direction of q. 
In many practical considerations one integrates over the azimuthal angles </> of the recoil nucleon, 
which yields: 



da 



x y m n 



dxdQ 2 d 3 p s /E s xQ 4 



Q 2 
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v 



F^ + 2tan 2 (-))- 
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(3) 



where: 



at s ,Pt) 



F L + 



2q 2 



mjv 



?D 
T 



and 



F^{x,Q 2 ,a S) p t ) 



The theoretical description of the reaction ([T]) at not very large values of p s < 700 MeV/c 
is based on the assumption that it proceeds through the interaction of virtual photon off one of 
the bound nucleons in the deuteron, while produced particles can interact in the final state with 
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(a) (b) 
Figure 1 : IA (a) and FSI (b) diagrams 



the other (spectator) nucleon. Since we are interested in x > 0.3 kinematics it is legitimate to 
neglect simultaneous interaction of 7* with two nucleons. Two main diagrams will contribute to 
the cross section of reaction ©: impulse approximation (IA) (Fig. la) and diagram representing a 
rescattering of the recoil nucleon off the products of DIS scattering (Fig. lb), which we will refer 
as final state interaction (FSI) diagram. 



2.1 Impulse Approximation (I A) 

In IA the recoil nucleon is a spectator of 7* scattering off the bound nucleon N. One can apply 
Feynman diagram rules to write down the IA amplitude in a formal form (see e.g. lfT2l[T3lO : 

Aj A =< X\J? m (Q 2 ,v,p s ) i d ~f s + m , 2 u{ V s)T d =< X\J^ m (Q\u,p / d -f s + m u(p s )T d 

(4) 

where T d is the covariant d — > pn transition vertex which is a smooth function at the pole of 
the struck nucleon propagator, and J^ m (Q 2 , v, x) represents the electromagnetic DIS operator of 
electron scattering off the bound nucleon. Here we suppressed the polarization indices of the 
deuteron and the nucleons. 

Taking the recoil nucleon on mass shell in Fig. 1(a) and using j> d — j> s + m ~ u(p d — 

spins 

Ps)u(p d — p s ) one can factorize the IA amplitude into two parts, consisting of the DIS current of 
the bound nucleon, «7j£ N =< X\J£ m (Q 2 , u,p s )u(p d — p s ) and the wave function of the deuteron. 
The latter is expressed through the Y d vertex function and the bound nucleon propagator. Since 
the deuteron wave function is not a Lorentz invariant quantity its determination depends on the 
reference frame in which the above factorization is performed. 

The factorization procedure in the Lab frame of the deuteron in nonrelativistic limit (corre- 
sponding to the equal time quantization) yields the nonrelativistic deuteron wave function which 
is the solution of the Schroedinger equation: 

«f a (p.) = ■ (5) 

2^(2vr) 3 £ s K - (Pd - Ps) 2 

It is normalized as / \^ d (p s )\ 2 d 3 p s = 1. This correspondence usually achieved in the virtual 
nucleon (VN) approximation in which the scattering is described in the LAB frame of the nucleus 
and electrons scatter off the virtual nucleon whose virtuality is defined by the kinematic parameters 
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of the spectator nucleon. In this case the form of the wave function is defined through the evaluation 
of the IA amplitude at the one-mass shell pole of the spectator nucleon propagator in the Lab frame. 
This yields an off- energy - shell state of the bound nucleon. 

In another approximation one can formally associate the T d vertex with Light Cone deuteron 
wave function by considering equal light cone time r = t — z quantization. In this case: 



LC, 



-M 2 



(6) 



normalized as / |\E> d(a s ,p s t)\ 2 d 2 p s t^r = 1. In this approximation referred as the light cone (LC) 
approximation ifTTI the scattering is described in the light cone reference frame where the wave 
function is evaluated at the pole of the spectator nucleon in the LC reference frame. This yields an 
off - light cone - energy (E + p z )-shell state of the bound nucleon. 

The above described factorization of the scattering amplitude into two parts in IA allows one to 
express the nuclear DIS structure functions through the convolution of bound nucleon DIS structure 
functions and the nuclear spectral function, S as follows: ll6l fT5l : 
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71 



F$(x,Q 2 ,a,p t ) + ^F*f/(x,Q 2 ,a, Pt ) 
2pq 



(7) 



where FfU and FoU are the structure functions of the bound nucleon and sin 2 5 



?eff 



IN 



2N 



%. The 
q 

nuclear spectral function, S describes the probability of finding an interacting nucleon in the target 
with momentum (a, p t ) and a recoil nucleon in the final state of the reaction with momentum (a s , 
p st ). Note that in IA a s + a = 2 and p t = — p s t- S is a model dependent quantity whose form 
depends on the framework/formalism used to describe the interaction. 

In VN approximation n = 2 (m^-e s ) anc * ^ or tne s P ectra l function one obtains lfT5ll : 



S^(a s ,p st ) = E s ^ 2 (p s 



(8) 



where / Sy J %{a s ^p st ^d 2 p t = 1. It can be shown fT51 that S VN has a normalization such that 
in Eq.© it conserves the number of the baryons involved in the reaction. On the other hand the 
momentum sum rule is not satisfied. This reflects the fact that due to the virtuality of the interacting 
nucleon part of its momentum is shared by the non-nucleonic degrees of freedom in nuclei. 
In LC approximation ifTTll n = (2 — a s ) and: 



^Lc( a siPst) 



^d C (a s ,p st ) 



Ek 



a. 



-*2(*). 



Eh 



Eh 



Pst 



(9) 



The remarkable property of the LC spectral function is that it fulfills both baryon number conser- 
vation (/ S^(a a ,Pt)^-cPpt = 1) and momentum sum rule (J aSj£(a a , p t ) ^d 2 p t = 1) (El) 
(for recent discussion see also IH31 ). 
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2.2 Final State Interaction 



Calculation beyond the impulse approximation is a very complex task in DIS. This includes FSI 
contribution in which a reinteraction is taken place between spectator nucleon and the products of 
DIS scattering off the struck nucleon. The typical FSI amplitude can be represented as in Fig. lb. 
Calculation of this requires a detailed understanding of the dynamics of DIS as well as its hadronic 
structure at the final state of the reaction. However the local character of DIS scattering off one of 
the nucleons in the deuteron allows us to represent FSI diagram of Fig. lb in the most general form 
as follows: 

d4 P s '/V„\A jem fr fl ^ fid ~ fis' + m N fa + 7Tl Ni 



A fsi = E / Jti^i^ s \ A ?si ■ G{X') • 



x , j i(2n) 4 . -l. (p d —pj) 2 -m% + iefi, -m 2 Nl +ie 

( 10 ) 

where J em (Q 2 , x) and Apsi represent operators of DIS and FSI scattering and G(X') is a notation 
for the propagation of intermediate state X' . The amplitude in Eq.dTUt is too general to allow any 
practical calculation. Our further discussion is based on the presence of at least one integration 
over the momentum of spectator nucleon in Eq.ifTUb. Complexity of the intermediate state X' and 
its final state interaction may lead to additional integrals over momenta in Eq.dTOl). Obviously 
for purpose of demonstrating the lack of singularity it is sufficient to consider the case of one 
integration over the momentum of the spectator nucleon. 



3 Analytic Properties at the Pole 

Next we discuss the analytic properties of IA (Eq.@) and FSI (Eq. dTOl) ') amplitudes at the pole of 
the struck nucleon propagator. 

3.1 I A Amplitude 

We first consider IA amplitude, in which the momentum of the struck nucleon is fixed by the 
kinematics of the spectator nucleon. It follows from Eq.@ that I A amplitude has a pole at t — 
m 2 N which corresponds to the struck-nucleon being on-mass shell. Using the following relation 
between t and kinetic energy of the recoil nucleon, T s : 

t = -2M d T s + m 2 N -\e d \(M d + m N -m s ) (11) 

it is straightforward to represent the amplitude of Eq.© in the following form: 

A IA =< X\J em (Q 2 ,x)\n > u(p d -p e )u(p a ) ^ \ j_ o M T ' (12) 

\e d \(M d + m N - m s ) + 2M d T s 

This equation shows that the pole is associated with the negative value of kinetic energy of the 
spectator nucleon, equal to the half of the magnitude of the deuteron binding energy e d . 

rppole ksj / ~i , m n ~ m P \ ^ \ e d\ ,-.~s 

2 1 + M d ~ 2 ' K } 
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Despite being in nonphysical region of the reaction dU) the very fact of the existence of the pole 
indicates that the extrapolation of the measured semiinclusive cross section to the positive region 
of t (negative kinetic energy region of the spectator nucleon) will allow us to isolate the "free" 
cross section of the electron - nucleon scattering. This can be true only if the singularity is unique 
to the I A amplitude. 



3.2 Loop Theorem and non IA Amplitudes 

Now we will prove that since all FSI type amplitudes contain at least one loop integral in the 
amplitude, they are regular at the one-shell pole of the struck nucleon propagator. 

It is sufficient to prove that no singularities exist for the amplitude that contains at least one 
loop involving the spectator nucleon, to prove that all higher order diagrams will have a smooth 
behavior at t — > m 2 N limit. In considering the FSI diagram, using Eq. lfTTT) . we will analyze the ana- 
lytic behavior of Apsi with respect to t. Due to the relation in Eq. (fTTT) all we need to demonstrate 
is that Apsi is a smooth function in T s — > — y limit. 

In Eq.([TTjl) we first evaluate the integral over d°p s > at the pole of the spectator nucleon: 
/ d P°s' J =^a +ie = —£jrj, then introducing k = p s - p s > one obtains: 

at N S 



A 



FSI 



d 3 k 



(X,s\A FSI G(X')r m (Q 2 ,x)\N)u(p d - Ps - k)u(p s - k) 
2(m N + T S - k Q ) 



x- 



-M 2 + 2M d (m + T S - ho) 



(14) 



To estimate the contribution of this amplitude at the pole of the struck nucleon propagator it is 
enough now to estimate the integral 



/ 



(15) 



M d \e b \+2M d T s -2M d h 



For relatively small k, using the relation: k = \Jm 2 N + p 2 s — ym 2 N + (p s — k) 2 w ^ — we 
see that in the limit T s — > — ^ integrand in Eq. (fT5t is finite, 



k 2 dk dk 



(16) 



This demonstrates that the integral in Eq. (fl4l is finite in T s — > — ^ limit . Note that the integral 
converges also in the large k limit since in this case the integrand is suppressed due to the ampli- 
tude of rescattering of the system X off the spectator nucleon. This result demonstrates that all 
contributions that contain at least one rescattering with spectator nucleon have a smooth behavior 
at the pole of the struck nucleon propagator as compared to the singular behavior of IA amplitude 
in Eq.©- 
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4 Extracting "Free" Neutron DIS Structure Function 



We demonstrate the procedure of extraction of "free" neutron structure functions using specific 
models for evaluation of the FSI. We then demonstrate that the result of the extrapolation is prac- 
tically insensitive to the parameters which define the strength of FSI in these models. The latter 
indicates the model independence of the procedure. 



4.1 Model for FSI 

We discuss a model in which the FSI of the spectator nucleon occurs coherently with the system 
produced in the deep inelastic electron- bound nucleon scattering (see also (16)). Such approach is 
rather well founded for energies relevant to Jefferson Lab[9|, in which the multiplicity of produced 
system is restricted to few hadrons. Due to the same reason a nonspectator contribution in which a 
nucleon with small momentum is produced in the 7* + N — > N + X process is strongly suppressed 
at x > 0.3 as it requires scattering off a nucleon with a very large momentum a > l + x. Therefore 
contributions to the cross section associated with the incoherent rescattering will not contain any 
terms singular at t — > m 2 N . Thus we expect that the main contribution to DIS cross section will 
come from the diagrams of type of Fig. lb in which FSI amplitude could be estimated in analogy 
to the processes of e + d—^e + p + n considered in Ref . lITSl [T31 fTSIl . In this approach one arrives 
at the similar eikonal form of the rescattering amplitude treating the effective cross section and the 
t-dependence of the rescattering as free parameters. In the calculation we neglect a small effect 
of modification of Fermi momentum of the struck nucleon in the 7*iV amplitude due to FSI. This 
is justified due to the fact that the dominant process of rescattering corresponds to the conserved 
light-cone fraction of the interacting nucleon lfT21[T51l (see also below). Within this approximation 
we can factorize 7*iV DIS amplitude from the FSI. This allows us to use the whole theoretical 
framework described in Sec l2.1l bv replacing S IA with the distorted spectral functions, S DWIA 
which include both IA and FSI effects. 

Using the procedure of the generalized eikonal approximation of Ref. lfT2l fT3ll within virtual 
nucleon approximation for the distorted spectral function one obtains: 



qDWIA 




d 2 kt 



f(kt)[* d (ps)-i*'(ps)] 



+ 



(17) 



where p s = (p sz ,Pst) 
W 2 = (q + M d - p s 



'sz 

) 2 and W$ 



{E s -m) + 



w 2 



Here 



(p sz - A VN ,p st - kt) and A VN = iMsLM t ^ 

(q + m) 2 . ^>' d results from the non-pole contribution in the 



integration over the longitudinal component of the transferred momentum (for details and for the 
expression of ty' d see Appendix B of Ref.[13|). Similar to the IA case an averaging over the 
polarizations of the deuteron is assumed. 
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In LC approximation [14] for distorted spectral function we obtain: 



oDWIA 




f(l t ) \¥ d (k)Mk) - i*l(k)K(k) 



f(h)[* d (k)-M'(k)} 



+ 



(18) 



wherefc = fe,fc t ) = (k z j^-A lc , k t -l t ) and A lc = M ^( ! ^)+ I1 i^. Here k represents the 
relative momentum of the target nucleons in the light cone reference frame. They are expressed 
through the light cone momentum fraction and the transverse momentum of the target nucleons 
according the relations given at the RHS of Eq.©. 

It follows from Eqs.dTTl) and (fT81 that in both cases the distorted spectral function has three 
distinctive terms. The first term represents the IA term discussed already in Sec. 12.11 the sec- 
ond term represents the interference between IA and FSI amplitudes which screens the 7*g? cross 
section, while the third term is the pure rescattering contribution. 

Note that one of the important differences between VN and LC approximations, is that at 
high energy limit (q, q ^> m), Avn is finite approaching to (E s — m), while Alc vanishes. 
This reflects the fact that at high energy limit the light cone momentum fraction of the interacting 
nucleons (such us a s ) conserves during final state interaction. 

The amplitude / in both Eqs. (fT7t and (fTSt describes the rescattering of the spectator nucleon 
off the products of deep inelastic scattering off struck nucleon. Since characteristic momentum 
transfer in FSI is relatively small (<g m) we can model this amplitude in the form of: f(k) = 
cr e ff(i + a)e~~ k ^ with cr e //, a and B being free parameters. For our numerical estimates we 
fix B = 8 GeV 2 and a = —0.2 corresponding to characteristic values of hadronic interaction at 
small momentum transfer. We will allow variation of cr e // between zero and 60-80 mb. The latter 
numbers represents that characteristic cross sections of scattering of the spectator nucleon off the 
"X" state containing one baryon and up to two mesons. 

The distorted spectral functions of Eqs.(fT7t and (fT8l now can be used in Eq.© to estimate the 
Fg and Ffp in the approximation which includes both IA and FSI effects. Such approximation 
we will refer as distorted wave impulse approximation (DWIA). 

In Fig. 2 we represent the predictions of VN and LC models for F^ calculated within PWIA 
and DWIA approximations. One can see from the figure that if the strength of XN interaction 
is varied within a reasonable range (cjxn < 60 mb) one finds a rather broad range of predictions. 
This highlights the necessity of the dedicated measurements |[l9j[20| of semiinclusive DIS reaction 
off the deuteron covering wide kinematic range of the recoil nucleon that will allow the focused 
comparison of the predictions of different theoretical approaches. It is important however that at 
a = 1 the difference between numerical results of VN and LC approaches is negligible leaving the 
source of the main uncertainty to the strength of the FSI. 



w 2 -w% 



9 




Figure 2: Angular dependence of for different values of recoil nucleon momenta. Dashed 
and solid curves correspond to PWIA and DWIA predictions within VN approximation. Curved 
labeled by circles corresponds to LC predictions. The total cross section of XN state in the FSI 
amplitude is taken as 60 mb. Arrows identify the angles corresponding to a = 1. 



4.2 Extraction Procedure for F2 

We discuss now the procedure that allows us to extract the "free" F 2 for the struck nucleon. One 
starts with the assumption that we don't know the strenght of the FSI for which we can assume 
various values of effective reinteraction cross section from to 80mb (see discussion above). 

To employ the extraction procedure based on the pole extrapolation method we introduce an 
extraction factor I(p s ,t) defined as follows: 

1 «-t) 2 1 



IiPs,t) 



E s [Res(V d (T pole ))Y 



pq 



1 + cos5) 2 {a + §>a q ) 2 + \sin 2 5 



(19) 



where Res{^ d{T po ie)) = ^ ^ 2mN GeV~^) where m N corresponds to the mass of the struck nu- 
cleon. Here C slightly depends on the particular potential used to calculate the deuteron wave 
function. For example C = 0.3939 for the wave function with Paris potential|21 ] and C = 0.3930 
for the Bonn potential! 22 1. 

Using I(p s , t) we define the extracted structure function as follows: 

F™ tr (Q 2 ,x,t) = I(p s ,t) ■ F% EXP (x,q 2 ,a s ,p t ), (20) 

SI EXP 

where F 2D ' is the experimentally measured value of the integrated structure function defined 
in Eq.®. 
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If only PWIA contributions were important in the cross section of the reaction ([T]) then one 
can estimate F2^' EXP (x,q 2 ,a s ,p t ) using Eq.© with spectral functions defined in Eqs.® and 
©. In this case one observes that by extrapolating t to the value of m N , F|^* r approaches to 
F^(x,Q 2 ,a = l,p t = 0) = F£T(x,Q 2 ). 

In the case when FSI effects are considered we evaluate F 2 s ^ ) ' EXP (x, q 2 , a s ,p t ) within DWIA 
framework, in which F^ EX p is estimated using Eq.© with the distorted spectral functions de- 
fined in Eqs. (fT7t and (fl8t . In this case we again extrapolate F 2D ' XP (x, q 2 , a s , pt) to the values of 
t -> m 2 N . 

To establish the appropriate extrapolation procedure one has to investigate the analytic proper- 
ties of Eq.(l2*0l with respect to the variable of t' = t — m 2 N . Here one observes that F|^ tr , at small 
\t'\ <C m 2 N , is an explicit quadratic function of t' with IA-FSI interference term being proportional 
to t', while the double scattering term (~ | Af,sj| 2 ) is proportional to t' 2 . 

Note however that the |A/a| 2 term has a (hidden) weak polynomial dependence on t' due to 
higher mass singularities in the deuteron wave function (ip d ~ f a with M 2 = (i— l)47m + 

i « 

2(i - l) 2 m° (see e.g. Ref.(23, where 7 = 0.04568 GeV and m = 0.17759 GeV. Since the second 
nearest pole corresponds to rather large positive values for t' (~ 0.1 Gev 2 ), based on Eq. (fTT1) on 
estimates that by restricting recoil nucleon momenta p s < 250 MeV/c the |Am| 2 term can be 
represented again as a quadratic function of t' . Therefore if spectator nucleon kinetic energy (or 
t') is much less than the energy scale corresponding to the higher mass singularity in the deuteron 
wave function, Eq. (l20l will represent a quadratic function of t' . 

The above discussion suggests that the pole extrapolation can be achieved by the square fit of 
the experimental data measured at small and finite values of t' . However such extrapolation can 
be meaningful only if we can exclude the variation of F|^ tr (Q 2 , x, a, t) due to the change of the 
other kinematic variables involved in the reaction. One significant variation is the a dependence 
of F extr due to the combination of - entering in the bound nucleon structure function. This may 
significantly alter t' dependence at large x kinematics. Another important effect in high x kinemat- 
ics is the higher twist effects due to sensitivity of structure functions on the produced final mass 
of the DIS scattering, Wn at intermediate Q 2 . Since W x ~ otWfj Q where Wno is me fi na l mass 
produced off the stationary nucleon, the Wn sensitivity will result in additional a dependence of 
F extr . Thus, in the large x kinematics, especially at intermediate Q 2 , the t'-dependence of F^ 
will change strongly with a. 

Since 

p s ■ cos(9 s ) 

a~l + — (21) 

m 

the discussed above sensitivity will result in the sensitivity of the extrapolation procedure to the 
direction of the momentum of the recoil nucleon. To illustrate this point we calculate the ratio, 



Fj N tr (Q 2 ,x,a, t) 
Fl r N ee (Q 2 ,x) 



= -27V (22) 



and study its variation with i! at different fixed values of 9 S . Figs. 3(a) and (b) demonstrate such cal- 
culations within both DWIA and PWIA using VN approximation at x = 0.7 for Q 2 = 5 Gev 2 (a) 
and Q 2 = 10 Gev 2 (b) kinematics. These calculations reveal significant violation of quadratic 
dependence of F extr at9 s > (<) 90° in which case the corresponding a is away from unity. This 
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Figure 3: The —(t — m\) dependence of R for different values of recoil nucleon angle. Dashed 
and solid curves correspond to PWIA and DWIA predictions within VN approximation. The total 
cross section of reinteraction of XN state in the FSI amplitude is taken 60 mb. 

reflects the fact that both F^ tr (Q 2 , x, a, t) ~ F^ nd (Q 2 , W N , t) and F( r N ee (Q 2 , x, W m ) are 
sensitive functions of x and W and a little mismatch between x and x w tt— as well as Wm and 
Wjvo can result to a substantial irregularity in t' dependence of R. Note that W sensitivity is espe- 
cially strong at small 6 S (see 6 S = 30° curves in Fig. 3a) which corresponds to small (resonating) 
values of Wn- However we can see from Fig. 3(b) that the W dependence diminishes with an 
increase of Q 2 due to suppression of the higher twist effects. 

Although both Fig. 3a and 3.b demonstrate substantial angular (a) dependencies, they also 
indicate the way to avoid these complications. Due to relation (ETT) these irregularities are gone at 
9 S = 90°, since in this case a ~ 1 and x ~ x and ~ W^ . Therefore one concludes that a s = 
a = 1 represents the most suitable kinematical condition for t' extrapolation of F|^ r (Q 2 , x, t). It 
is worth mentioning that according to Fig. 2, a s — 1 also diminishes the discrepancy between VN 
and LC approaches. 

Fig. 4 represents one example of the extrapolating procedure for a s — 1 kinematics. In this case 
R is calculated using both PWIA and DWIA models within VN approach for x = 0.8 and x = 0.9 
at Q 2 = 10 GeV 2 and for physically accessible values of t' < 0. Then we use three calculated 
points of R to make a quadratic fit which is then extrapolated to the t' — > values. For DWIA 
the FSI strength is set to 80 mb which may be considered somewhat above the highest possible 
value for the cross section of the X-state scattering off the spectator nucleon. For the fitting points 
we choose three R values calculated at — t' (0.1, 0.2, 0.3) corresponding to (54MeV/c, 89MeV/c, 
1 14MeV/c) values of recoil proton momenta. The choice of these three values is motivated by the 
kinematic acceptance of the experiment of Ref . ||9l . 
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Q 2 =10 GeV 2 , x=0.8 g Q 2 =10 GeV 2 , x=0.9 




-(t-m 2 ), GeV 2 -(t-m 2 ), GeV 2 

Figure 4: Extrapolation of R at — (t — m^) — > based on the square fit using first three points of 
calculated R. The triangles and open circles correspond to R (Eq. (l22l) calculated within PWIA and 
DWIA using VN model. 

The quadratic function fitted to above points yields the following extrapolated values for R at 
-if = 0: (0.9923 and 0.9888) for (PWIA and DWIA) predictions at x = 0.8 and (0.9868, 0.9828) 
for x = 0.9. These numbers demonstrate that extrapolation procedure practically eliminates the 
uncertainty due to the FSI, making FSI effects less than 0.5%. 

Similar estimations within LC approximation yield for extrapolated values of R at —if = 0: 
(0.9983, 0.9946) for (PWIA, DWIA) predictions at x = 0.8 and (1.0078,1.0036) for x = 0.9. 
These numbers combined with the numbers obtained above in the VN approximation demonstrate 
that the overall uncertainty of extrapolation procedure using the fitting points consistent with the 
kinematics of experiment of Ref.[9| is of the order of 1%. This is an unprecedented accuracy that 
can be achieved in extraction of high x DIS structure function of the nucleons. 

Note that due to the higher order singularities in the deuteron wave function, moving the fitting 
points away from the t' — limit will worsen the convergence of the quadratically fitted function 
of R to the unity. In this case an additional improvement of the procedure can be achieved if we 
additionally normalize R in Eq.(l2"2l by the momentum distribution of the deuteron. 

5 Conclusions 

We have considered the pole extrapolation procedure aimed at model independent extraction of 
high x DIS structure function of struck nucleon in semiinclusive DIS scattering off the deuteron. 
For analysis of the validity of this procedure we considered two theoretical approaches (virtual 
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nucleon and light cone) for description of semiinclusive deep inelastic scattering from the deuteron. 
Using these approaches we modeled the final state interaction based on distorted wave impulse 
approximation. 

Within this framework we demonstrated that a simple quadratic fit is well suited for on-mass 
shell extrapolation of the bound nucleon structure function. Although the strength of FSI is varied 
in the wide range of the scattering cross sections 4- 80 mb the overall uncertainty due to FSI 
and the choice of the particular (VN or LC) theoretical approximation is estimated to be on the 
level of 1%. This result indicates the model independent character of the mass-shell extrapolation 
procedure. 

The numerical estimates are done for t' values characteristic to the experiment of Ref.[9|. 
Thus one can expect the same level of uncertainty in the extraction of on-shell structure functions 
of nucleon if the pole extrapolation procedure is applied to the actual data. 
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